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Abstract 

In this paper, we consider the second-order equations of Duffing 
type. Bounds for the derivative of the restoring force are given that 
ensure the existence and uniqueness of a periodic solution. Further- 
more, the stability of the unique periodic solution is analyzed; the 
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sharp rate of exponential decay is determined for a solution that is 
near to the unique periodic solution. 
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AMS subject classification: 34B15, 34D05. 

1 Introduction and statement of main results 

This paper is devoted to the existence, uniqueness and stability of periodic 
solutions of the Duffmg-type equation 

(1.1) x" + cx' + g(t, x) = h(t), 

where 3: R x 1 ^ M is a T-periodic function in t and continuous in x, 
h{t) is a T-periodic function, and c > is a positive constant. The exis- 
tence and multiplicity of periodic solutions of (jl.ip or more general types 
of nonlinear second-order differential equations have been investigated ex- 
tensively by many authors. However, the stability of periodic solutions is 
less extensively studied. In [9], R. Ortega studied (jl.ip from the stability 
point of view and obtained stability results by topological index. A.C. Lazer 
and P.J. McKenna get stability results by converting the equation (jl.ip to 
a fixed-point problem [5]. Recently, more complete results concerning the 
stability of periodic solutions of (jl.ip were obtained by J.M. Alonso and 
R. Ortega [TJ [2] . Under the condition that the derivative of the restoring 
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force is independent of t and positive, they found sharp bounds that guar- 
antee global asymptotic stability. In [2], optimal bounds for stability are 
obtained. However, the rate of exponential decay of periodic solutions is 
not known in general. The purpose of this paper is to determine the rate of 
decay of the periodic solutions. Our results are motivated by careful obser- 
vation of the following example of a second-order differential equation with 
constant coefficients: 



where c > 0. If k ^ 0, then (jl.2p has a unique T-periodic solution xo(t), 
and all solutions of (|1.2|) have the form x(t) = c\e plt + C2e P2t + xo(t) where 
pi < p2 are roots of p 2 + cp + k = 0. Regarding A; as a parameter, when 
k < 0, then p2 > 0, and the unique periodic solution is unstable; when k > 0, 
the unique periodic solution is stable; keep increasing k till it crosses the 
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critical value k = -r , then p\, p2 are a pair of conjugates. In this case, every 
solution other than the unique T-periodic solution decays to the periodic 
solution at the same exponential rate ^ independently of k. 

The aim of this paper is to show that the above-mentioned | decay results 
can be generalized to (II. ip . In order to show that periodic solutions of (jl.ip 
are stable, it is essential to impose conditions on the restoring force that can 
rule out the existence of additional periodic solutions that are subharmonic 
of order 2, whereas to determine the rate of decay of periodic solutions, 




x" + cx' + kx = h(t) 
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it is essential to impose conditions on g(x,t) such that the corresponding 
Dirichlet boundary- value problem does not admit any nontrivial solutions. 
We shall give existence and uniqueness results, characterizing solutions that 
are locally asymptotically stable with sharp exponential rate of decay, by a 
Sturm comparison argument coupled with Floquet theory. Roughly speak- 
ing, (jl.ip has a unique T-periodic solution which is locally asymptotically 
stable in the sense of Lyapunov, as long as the restoring force is small or the 
fractional constant is large. 

The following notations will be used throughout the rest of the paper. 

1. lP T T-periodic function u € T p [0, T] with ||w|| p for 1 < p < oo; 

2. T-periodic function u G C k [0, T], k > 0, with C k -norm; 

3. a(t) ^> P(t), if a(t) > (3(t) and a(t) > (3(t) on some positive-measure 
subset. 

Now we state our main results. 

Theorem 1.1. Assume that g(t,x) G C 1 (1R) satisfies the following condi- 
tions: 

(1) g' x (t,x) <^z + ^, and 
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(2) there is an a(t) € Ct with a(t) 3> such that g' x (t,x) S> a(t), where 
a(t) denotes the the average of a(t) over a period. 
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Then (jl.ip has a unique T -"periodic solution which is asymptotically sta- 
ble with sharp rate of decay of | for c > 0. 

Here, we say that the periodic solution uq of (11. lj) is locally asymptot- 
ically stable if there exist constants C > and a > such that if -u is 
another solution with \u(0) — uo(0)\ + \u'(0) — u' Q (0)\ = d sufficiently small, 
then \u(t) — uo(t) | + \u'(t) — u' Q (t) \ < Cde~ at . The super exponent a as above 
is called the rate of decay of uq. 

Remark. The bounds given in Theorem ll.ll are optimal. In fact, the example 
given at the beginning of the section shows that the lower bound is sharp. 
An example to show that the upper bound is optimal in the theorem will be 
given in the following. 

Example. Consider the linear differential equation 

(1.3) x" + cx' + ^(l + c 2 + ecost)x = 0, 

where c > and |e| is small and e / 0. By the transformation y(t) = 
e~2*x(t), the damping term can be eliminated and equation (II. 3D is equiv- 
alent to 

(1.4) y" + ^l + ecost )y = 0. 

The following notations are used. Obviously, the nontrivial solution x(t) of 
(jl.3p has rate of decay | if and only if the corresponding solution y(t) of 
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(jl.4p is a nontrivial bounded solution. Let y^, i = l,2, denote the solutions 
of (|1.4p with initial values 

Vl (0, e) = y' 2 (0, e) = 1, ^(0, e) = y 2 (0, e) = 0. 

If the discriminant function of (|1.3|) is denoted by A(e) = yi(27r, e) + 
2/2 (2?r, e), then the Floquet multipliers are the roots of the quadratic 

// 2 - A(e)/i + 1 = 0. 

By means of perturbation one can compute that 

A( e ) = -2-^e 2 + 0(e 2 ) 

for e small, which shows that the modulus of one of the multipliers is greater 
than 1 and the modulus of the other multiplier is less than 1. Therefore the 
multipliers of (|l-4[) are a pair of distinct real numbers. This implies that 
the rate of decay of (|1.4p is greater than |. Thus if the derivative of the 
restoring term crosses over the bound given in Theorem II. 1\ the conclusion 
of Theorem 11.11 does not hold. 

Consider the piecewise linear equation 

(1.5) x" + cx' + a{t)x + - b{t)x~ = h(t), 

where x + = max{x,0}, x~ = min{— x,0}, a(t),b(t) G Ct- For a and b 
constant, (|1.5|) is a simple mathematical model for vertical oscillations of 
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a long-span suspension bridge, which has received great attention after a 
series of works of Lazer and McKenna [6] . 

Theorem 1.2. Suppose that h(t) E Ct has a finite number of zeros in 
[0,T]. Then 



(1) equation (|1.5p has a unique T-periodic solution if 
0«a(t),fe(t)«^ + $, 

(2) the unique T-periodic solution is asymptotically stable if 
0<a(i),6(i) + and 

(3) the unique T-periodic solution has rate of decay % if 
£ <a(t),&(*)«&£ + T- 



2 The linear periodic problem 

In this section we shall recall some basic results about topological methods 
and prove some stability results for linear periodic systems. 
Consider the periodic boundary-value problem 



(2.1) 



x' = F(t,x), 
x(0) = x{T). 

where F : [ 0, T ] x M. n — > M. n is a continuous function that is T-periodic in 
t. In order to use a nomotopic method to compute the degree, we assume 
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that h: [0,T] x W l x [0, 1] — >• R n is a continuous function such that 

h(t,x,l) = F(t,x), 
h(t,x,0) = G(x), 

where G(x) is continuous. The following theorem is due to J. Mawhin [8]. 

Lemma 2.1. Let £1 C Ct be an open bounded set such that the following 
conditions are satisfied. 

(1) There is no x E dQ such that 

x = h(t, x, A) VAG[0,1). 

(2) deg{g,nnR n ,0) + 0. 

Then (|2.ip /ias at least one solution. 

Next we consider the above system for n = 2. We denote by x(t,xo) 
the initial- value solution of (|2.ip and introduce the Poincare map P: xo — > 
x(T, xo). It is well known that xo) is a T-periodic solution of the system 
(|2.ip if and only if xo is a fixed point of P. If x is an isolated T-periodic 
solution of (|2.ip . then xo is an isolated fixed point of P. Hence the Brouwer 
index is defined by 

ind[P,x ] =deg(/-P,P e (0),0). 

Definition. A T-periodic solution x of (|2.ip will be called a nondegenerate 
T-periodic solution if the linearized equation 

(2-2) y' = f x (t,x)y 
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does not admit any nontrivial T-periodic solutions. 

Let M(t) be the fundamental matrix of (|2.2|) and fi\ and \ii the eigen- 
values of the matrix M(T). Then x(t, xo) is asymptotically stable if and 
only if < 1, i = 1, 2: otherwise, if there is an eigenvalue of M(T) with 
modulus greater than one, then x(t, xo) is unstable. 

Before giving some results concerning stability of linear periodic equa- 
tions, we consider the following eigenvalue problem: 

x" + Ax = 0, 

x(0) = x(T) = 0, sgnx'(O) = sgnx'(T). 

It is easy to verify that A n = ^ 2r ^ , the n-th eigenvalue of (|2.3p with the 
eigenfunction ip n (t) = sin^S^t. Now for 



(2.3) 



(2.4) 



x" + q(t)x = 0, 

x(0) = x(T) = 0, sgnx'(O) = sgnx'(T), 
we have the following result. 



Lemma 2.2. Assume that q(t) <C ^} . Then equation (|2.4p does not admit 
any nontrivial solutions. 

Proof. Suppose that x(t) is a nontrivial solution of (|2.4p . Then x is an 
eigenfunction. There is an n such that X n (q) = 0. Since q(t) <C , 
by a theorem concerning comparison of eigenvalues, we have that A n (g) > 
A n = A n — Ai > 0, a contradiction. □ 



Consider the homogeneous periodic equation 
(2.5) L a x := x" + cx' + a(t)x = 0, 

where c £ K is constant and £ L^f . 

Lemma 2.3. Assume a(t) £ L^? satisfies the following conditions: a(t) -C 
£p£ + ^andMt)>0. 

Then (|2.5|) does no£ admit any nontrivial T '-periodic solutions. 

Proof. Suppose on the contrary that (|2.5p admits a nontrivial T-periodic 
solution We claim that x(t) vanishes at some to £ [0,T]. If not, 

then x(i) 7^ for all t in R. By the periodic boundary conditions, we have 
x'(T) = x'(0) and x Jfr} = Dividing (|2.5p by and integrating by 

parts gives that 

f T x'(t) 2 f T 
J lW dt + Jo a{t)dt = °> 

which contradicts the hypothesis of the lemma. So x(t) has a zero in [0, T]. 
We may assume that x(0) = so that x(0) = x(T) = 0. By the transforma- 
tion y(t) = e2*x(t), y(t) is a nontrivial solution of (|2.5p with q(t) = a{t) -j, 
and if the first condition of the lemma holds, q(t) <C ^3 . Then according 
to Lemma 12.21 y(t) = 0, hence x(t) = 0, a contradiction. □ 

The following simple lemma, given in [4], will be used in proving the 
existence and the uniqueness of periodic solutions. 
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Lemma 2.4. Suppose that a{t), a\(t) and 0:2 (i) € L^? such that a.\(t), 



(1) f/ie possible T -periodic solution x of equation ()2,5p is either trivial or 
different from zero for each t £ R; 

(2) L ai x = (i = 1,2) cannot admit nontrivial T -periodic solutions si- 
multaneously if ot\{t) -C 02 (i); and 

(3) L a x = /ias no nontrivial T-periodic solution, if a(t) 3> (resp., 
a(t) <C 0). 

The following lemma is essential for determining the rates of decay of 
periodic solutions of (jl.ip . 



Lemma 2.5. Assume that a(t) <C ^ + ^ and i/iai a(t) > T/ien 
does no£ admit real Floquet multipliers. 

Proof. If the conclusion of the lemma does not hold, then there is a real 
Floquet multiplier p and a nontrivial solution x(t) such that x(i+T) = px{i). 
Introduce the transformation x = e~z ct u. Then u solves the equation 



(2.6) u" + 



«(*) " J 



u = 



with Floquet multiplier x = e 2 P- If x does not change sign in R, neither 
does u. Dividing equation f)2.6|) by n and integrating by parts, and noting 
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that u j£} = ^-t^t, we have that 



u(T) u(0) 



o u{t) 2 



<*(<) - j 



dt = 0, 



which contradicts the assumption that a(t) > Therefore x(t) vanishes at 
some to ^ [Oi^! an d the conditions are such that x(to + T) = px(to) = 0. 
Thus the corresponding u is a nontrivial solution of the following Dirichlet 
boundary-value problem: 



(2.7) 



<*(t) - — 



u = 0,u(t ) = u(t + T) = 0. 



2 2 

Since a(t) — <C |?r, multiplying (|2.7p by u and integrating we have that 



T+t 



u' 2 dt -- 
to J to 



T+to 



<*(t) - — 



IT 2 f T+t0 

u 2 dt < — g / u 2 dt, 
T J ^ 



which contradicts the Poincare inequality. Therefore (|2.5jl does not admit 
real multipliers. □ 



Lemma 2.6. Under the conditions of Lemma I2.5L the rate of decay of any 



nontrivial solution of f|2.5|) is S. 



Proof. Consider the corresponding system 



(2i 



X'(t) = A(t)X(t), 



where the column vector function X(t) = (x(t),x'(t)) T and -A(t) is the ma- 
trix function 



A(t) 



I o i N 



\ ~P(t) ~C J 
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Let M(t) be a fundamental matrix solution of (|2.8p . It is well known that 
M(t) has the form 

(2.9) M(t) = P(t)e Bt , 

where P(t) and B are 2x2 matrices, P(t) = P(t + T), and B is a constant 
matrix. Let p\ = e TXl , p2 = e TX2 be the Floquet multipliers and Ai and 
A2 the Floquet exponents associated with p\ and p2- Let x\ and X2 be the 
eigenvector components of the matrix e TB . It follows from Lemma 12.51 that 
pi and P2 are a pair of conjugates. Thus the eigenvectors associated with 
different eigenvalues are linearly independent. Therefore yi = pi(t)e Xit (for 
i = 1,2) form the fundamental solution system of equation (|2.8p . On the 
other hand, by applying the Jacobi-Liouville formula, we have 

\ pi f = P1P2 = e - fo T ^ = 

and 

Re Ai = Re A 2 = - Re (Ai + A 2 ) 

Since every solution is a linear combination of yi(t) and y2(t), Pi(t) is T- 
periodic, hence it is bounded. Therefore every nonzero solution of equation 
(I2.5h decays at the same exponential rate of 5. □ 
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3 Proof of main results 

Now we prove our main results. 
3.1 Proof of Theorem fTTTI 

We begin with the following existence result, dividing the proof into two 
steps. 

Step 1. Existence. 
Define F : Cf, -» C T by 

F(x(t)) := x" + cx' + g(t, x(t)). 

We have the following. 

Theorem 3.1. Assume that g(x,t) E C 1 (M x M), and that g(x,t) is T- 
periodic in t. If in addition g satisfies the two conditions 

(1) g(t,x)/x < + f Vx el and 

(2) there is a T -periodic function (3{t) G Ct such that (5{t) > and 
g(t,x)/x » j3(t) for all x G R, 

then the differential equation (II. ID has a T -periodic solution. 

Proof. Without loss of generality, we may assume that g(0, t) = 0, for oth- 
erwise we can reduce both sides of equation (jl.ip by g(0,t). Consider the 
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parametrized equation 

(3.1) Fx := x" + cx' + Xg(t, x) + (1 - \)ax = Xh(t) 



for some a G ^0, — h ^ ) . We claim that there is an R > such that 
equation (|3.ip has no solution on for any A € [0,1]. If there is not 
such an R, let x n be a sequence of solutions such that \\x n \\ — > oo and 
A n G [0,1]. Denote by z„ the quotient z n := jjf 12 ^- Dividing (|3.ip by ||x n ||, 
then multiplying by tp(t) G and integrating by parts, we have that 

(3.2) / z n ip" - cz n ip' + g(t,x n )ip/\\x n \\ dt = X n / iph n /\\x n \\ dt. 

Jo Jo 

The conditions of Theorem l3.1l implv that {[\ n g(t, x n )+(l—X n )a]x n /\\x n \\} 

is bounded. It is pre-compact in the weak topology in L 1 [0,T]. Thus there 

are subsequences such that g(t, x n )/x n — > g(t) and A n — > A. Taking the limit 

in the equation (|3.2p . one obtains that 



2 c 2 



(3.3) / {z^" - czip' + w(t)<pz} dt = 0, 

Jo 



where u;(t) = Xq(t) + (1 - A)a. Evidently, > Aa(t) + (1 - A)a > 0, 

and w(t) < — h ^, which satisfies the conditions of Lemma 2.2 in [3], 
and since here z(t) is a T-periodic solution, it follows from Lemma 12.31 that 
z(t) = 0, which contradicts the fact that, by construction, ||.z(i)| = 1. Next, 
by applying the homotopic invariance property, we have that 

deg(F ,B R ,0) = deg(F u B R ,0) = 1. 
15 



According to Lemma [2. 1\ equation (jl.ip then has a T-periodic solution. This 
completes the proof of existence. 

Now, we will finish the proof of Theorem 13.11 by verifying uniqueness. 
Let x\ and x^ be two distinct T-periodic solutions of the equation (jl.ip . and 
let u = x\ — X2- Then u satisfies the equation 



where p(t) = [g(t,xi(t)) — g(t, xi(t))]/[x\ — x%[. If u ^ identically, it will 
imply that u is an eigenfunction associated with a Floquet multiplier equal 
to one. Again, Lemma 12.31 rules out this possibility. Therefore u = 0. □ 

Step 2. Rate of decay. 

To show that every solution of the nonlinear equation (jl.ip locally decays 
at the rate of | to the unique T-periodic solution, we need the following C 1 
version of the Hartman-Grobman theorem [3] . 

Lemma 3.1. Let f : U C R n -> R n be a C 1 function with /(0) = such 
that fx(fi): R n — > M n is a contracting mapping. Then f is C 1 conjugate 
equivalent to /^(0). 

Consider the planar system associated with equation (11. ip . 



(3.4) 



u" + cv! + p(t)u = 0, 



x = cx — y, 



(3.5) 



< 



y 



h{t)-g(t,x). 
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Let Xo(t) = (xo(t),yo(t)) be the unique T-periodic solution determined 
by the initial condition Xo(0) = (xo> 2/o)- Then Xq corresponds to the unique 
fixed point of the Poincare mapping PX = U(T,X), where U(t,X) is the 
initial-value solution of (|3.5p with U(0, X) = X. Let M(t) be the funda- 
mental matrix solution of the linearization 

(3.6) X' = A(t)X 



of (|3.5p . where 

I a -1 

A(t) = 

\ -Q(t) 

By the differentiability of X(t) with respect to the initial value, the Poincare 
mapping can be expressed in terms of the initial value X by the following 
formula: 

(3.7) PX-X = M(T)(X-X ) + o(X-X ). 

Referring to Lemma 12.51 M(T) has a pair of conjugate eigenvalues A, 
A with |A| = e" cT / 2 . Thus P{X) is a contracting mapping. According 
to Lemma 13. 1\ there is a C 1 diffeomorphism (p which is near enough to 
the identity that PX — Xq is conjugate equivalent to M(T). There is an 
invertible constant matrix C such that 



C- L M{T)C 



A 



D(X), 
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and we may suppose that 

(3.8) i \X - X \ < \tp(X) - <p{X )\ <2\X-X \ 

for X — Xq small, since (p is near the identity. Therefore, the Liapunov 
exponent is given by 

ii x = lim -Llnli^X-Xol 

ra— >oo nl 

= lim — In 1^ o M(T) n o y~ l {X) -ipo M(T) n o y _1 (Xo)| 
= lim -llnl^AfoC- 1 [^(X)-^- 1 ^)]^-^. 

n^oo nl '2 

The second equality follows from (|3.8p . Hence, the rate of decay of the 
solution to the unique T-periodic solution is c/2, independently of the initial 
value X. 

Remark. The above result shows that the Lyapunov exponent is invariant 
under a C 1 conjugate transformation. From the proof, the conclusion is still 
true if ip is a bi-Lipschitz mapping. 

3.2 Proof of Theorem fTT2l 

The existence and uniqueness of T-periodic solutions of (|1.5j) can be obtained 
by the same argument in the proof of Theorem 11.11 But for the stability 
of the periodic solution, the C 1 regularity of the solutions of (|1.5H with 
respect to the initial value is needed. The following result is due to Lazer 
and McKenna. 
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Lemma 3.2. Let u(t, £, rj) be the solution of the initial-value problem 

u" + cv! + a(t)u + - b{t)u~ = h(t), 
^ «(0) = £, u'(Q) = n. 
Assume that a(t), b(t) G Ct, and that h G Ct has a finite number of zeros. 
Then, for t G [0,i], the partial derivatives of u and u' with respect to £,rj 
exist and are continuous. Moreover, if 



X(t) 



tin 


du 


ae 


dr) 


du' 


du' 


dt 


drj 



then 



X'(t) = A(t)X(t), X(0) = Id, 



A(t) 



where 

( 

1 

-pit) -c 

and p[t) = a(t)x+ +b(t)x-, where x± denotes the characteristic function of 
the set {t G [0,T], u(t) ^ 0}. 

Remark. Though Lemma 13.21 was originally stated only for both a and b 
constant, by carefully examining the proof in the appendix of [7], one may 
assert that Lemma 13.21 holds when a(t) and b(t) are continuous. 

Now that u(t) and u'(t) are known to be C 1 with respect to the initial 
value, by reasoning along the same lines as in the proof of Theorem II .1| the 
stability and rate of decay of the periodic solutions can be obtained. 
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